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Complex-space singularities of 2D Euler flow in Lagrangian coordinates 
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We show that, for two-dimensional space-periodic incompressible flow, the solution can be evalu- 
ated numerically in Lagrangian coordinates with the same accuracy achieved in standard Eulerian 
spectral methods. This allows the determination of complex-space Lagrangian singularities. La- 
grangian singularities are found to be closer to the real domain than Eulerian singularities and seem 
to correspond to fluid particles which escape to (complex) infinity by the current time. Various 
mathematical conjectures regarding Eulerian/Lagrangian singularities are presented. 
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I. INTRODUCTION 

Solutions to the incompressible Euler equation, start- 
ing from entire initial data (e.g. trigonometric polynomi- 
als), can be analytically continued to the complex space 
as long as they stay analytic in the real space. Further- 
more it is known since the seventies that any singularities 
in the real space, if they exist.have to be preceded by 
complex-space singularities 0, | . In 2D Euler flows, it 
is known that initial real-space analyticity for periodic 
solutions is never lost. This was proven in Refs. [1,01 hi 
which it was shown that the distance 8(t) to the real do- 
main of the nearest complex-space singularity, measured 
by the exponential falloff of the Fourier amplitude, de- 
creases at large times at most as a double exponential. 
Actually, already twenty five years ago, spectral simula- 
tions with 256 2 Fourier modes indicated that the decrease 
is more like a simple exponential Q . Spectral simulations 
at much higher resolutions, up to 8192 2 , which will be 
published elsewhere also indicate a behavior much closer 
to a single than to a double exponential. 1 The discrep- 
ancy between the simple and the double exponential is 
generally believed to be due to the phenomenon of de- 
pletion: the flow organizes itself into ribbon- like vortical 
structures in which the nonlinearity is almost completely 
suppressed (the nonlinearity would vanish identically if 
the flow depended on a single Cartesian coordinate). The 
same phenomenon also exists in three dimensions and 
could conceivably prevent finite-time blow up. 

In three dimensions the Beale-Kato-Majda (BKM) 
theorem implies that any blow up must be accompanied 
by the unboundedness of the modulus of the vorticity 
in the real domain Q (see also Q). In two dimensions, 
when the initial vorticity is bounded, this is of course 
ruled out by vorticity conservation. More precisely, it is 
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1 We nevertheless conjecture that by suitable analytic regulariza- 

tion of the flow considered in Ref. [g long-lasting transients with 

a double exponential decrease can be achieved. 



ruled out in the real domain, but not in the complex do- 
main. Actually, increasingly strong numerical evidence 
has been obtained indicating that the vorticity is infi- 
nite at complex singularities [iiGilEU- 2 Such numerical 
results were obtained only for flows in which the initial 
stream function is a trigonometric polynomial (the 2D 
analogues of the famous Taylor-Green flow [H|]), which 
are instances of entire functions, that is, analytic func- 
tions that have no singularity at finite complex locations. 

The "experimental result" about infinite vorticity 
along the complex singularities in two dimensions has 
an important consequence: because the conservation of 
vorticity along fluid particle trajectories carries over to 
complex trajectories, (Eulerian) complex locations with 
infinite vorticity are associated with fluid particles ini- 
tially at complex infinity; indeed, this is the only place 
where an entire function can be infinite. We were thus 
led to investigate the issue of (complex) singularities in 
Lagrangian coordinates. A Lagrangian singularity is a 
location at which the (analytic continuation of the) La- 
grangian map goes singular. Could it be that for two- 
dimensional flow there are no (complex) Lagrangian sin- 
gularities at finite distance? In other words: does the 
flow in Lagrangian coordinates preserve its initial entire 
character? A few years ago we performed very accurate 
numerical simulations, reported here for the first time, 
and we found strong evidence that the answer is "no" . 
W. Pauls and one of us (TM) [l4| then found a very 
simple counterexample to the preservation of the entire 
character: the "AB flow" ip = sinjricosa^ is an entire 
steady solution to the 2D Euler equation in Eulerian co- 
ordinates. For this flow, the trajectories of fluid particles 
can be expressed by elliptic functions and it was shown 
that, for any real positive time t, there exist complex 
initial locations of fluid particles which are mapped to 
infinity at time t and which thus are Lagrangian singu- 
larities. 



2 We have tried but failed to derive such a result from a complex 
version of the the BKM argument. 
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There is a considerable renewal of interest in the La- 
grangian structure of flows, both from a theoretical and 
experimental point of view (such issues frequently came 
up during the Euler conference). It is thus of interest to 
show that the Lagrangian description of flows can be ob- 
tained numerically with an accuracy comparable to that 
available by spectral methods for the Eulerian descrip- 
tion. The present paper is organized as follows. In Sec- 
tion [IT] we describe two numerical algorithms, which can 
be used for Lagrangian integration. In Section HTll we ap- 
ply this to the identification of complex Lagrangian sin- 
gularities. Here, all numerical studies are presented for 
the (unsteady) 2D flow with the simple initial condition 

ijjQ = COSXi + COS 2X2 , (1) 

which has been used in Refs. i, El EDI; key results arc 
also checked with the flow 

ipo = cos xi + cos 2x2 + sin(2a;i — X2) , (2) 

which has less symmetry than (JXJ) . Some concluding re- 
marks, with emphasis on mathematical conjectures, are 
presented in Section ITVl 



II. NUMERICAL SOLUTION IN LAGRANGIAN 
COORDINATES WITH SPECTRAL ACCURACY 

Our goal here is to obtain the velocity field u as a 
function of the Lagrangian location a and the time t. 
This Lagrangian field will be denoted u L (a, t). 

With simple boundary conditions, e.g. spatial period- 
icity, the easiest way to obtain high accuracy in a Eu- 
lerian simulation is to use a spectral or pseudo-spectral 
method [l3j . For analytic flow, whose Fourier transform 
decreases exponentially at high wavenumbers, the trun- 
cation error will then also decrease exponentially with 
the resolution. 

How does one carry this over to Lagrangian coordi- 
nates? In principle one can write an integro-differential 
equation for the (time-dependent) Lagrangian map a 1— > 
x. This equation has however nonlinearities with denom- 
inators which are not easily handled numerically. 

We present here two alternative methods, the spectral 
particle-tracking method (Section III A[) and the spectral 
displacement-Newton method (Section III B|) . 



A. Particle tracking method 

Obviously, the Lagrangian velocity field can be ob- 
tained by composing the Eulerian velocity field u(x, t) 
with the Lagrangian map x(a, t). The former can be ob- 
tained by standard spectral integration. The latter is the 
solution of the characteristic equation 

dtx(a, t) = u(x(a, t), t), x(a, 0) = a . (3) 



In the tracking method, we select a uniform grid of La- 
grangian points and "track" the fluid particles by inte- 
grating ((3]) along all the relevant fluid particle trajecto- 
ries. This can be done, e.g. using a fourth-order Runge- 
Kutta method. The problem is that, even if the initial 
positions coincide with Eulerian collocation points, this 
usually ceases to hold subsequently. Hence the Eulerian 
field must be interpolated. In order not to loose the spec- 
tral accuracy, the interpolation can be done using the 
Fourier series representation 

u{x, t) = ^w(fc,t)e ifc x . (4) 

k 

A difficulty is that, since the relevant x's are not collo- 
cation points, the velocities given by $4$ cannot be eval- 
uated using fast Fourier transforms but must be calcu- 
lated "naively" in 0(N 4 ) operations if we use an N x N 
grid. Furthermore this has to be done at every time step. 
Since the number of time steps needed to reach a given 
time t order unity is proportional to the resolution N, 
this method has a fairly large computational complexity 
0(N 5 ) and thus also a significant accumulation of round- 
off errors. For large values of the resolution N (512 or 
more) the particle tracking method is not very practi- 
cal unless we restrict the Lagrangian grid to being much 
coarser than the Eulerian grid. 

B. Displacement-Newton method 

This method makes use of the fact that the inverse 
Lagrangian map a(x, t) satisfies, in Eulerian coordinates, 
the equation 

d t a + u{x,t) ■ Va = , (5) 

which just expresses the constancy of the Lagrangian lo- 
cation a under advection by the velocity field. This equa- 
tion can be solved along with the basic Euler equation, 
both in Eulerian coordinates. This will however yield a 
map which still has to be inverted to obtain the direct 
Lagrangian map. 

For periodic boundary conditions the direct and in- 
verse maps are not periodic and it is more convenient to 
work with the displacement field, here defined as 

d(x, t) = a(x, t) — x . (6) 

It follows from ([5]) and ^ that the displacement satisfies 
the following equation in the Eulerian coordinates 

d t d(x, t) + (u(x, t) • V)d(jc, t) = -u(x, i), (7) 

with the initial condition d(x,0) — 0. This equation can 
be solved along with the Euler equation to obtain the 
displacement in Eulerian coordinates on a uniform grid 
of N x N collocation points. 

Then comes the difficult step, namely the inversion. 
For this we define the off-grid displacement, as above, by 
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its Fourier series, extended off-grid and we try to find the 
x locations associated to a set of Lagrangian collocation 
points on the regular grid A = (2iri/N, 2nj/N), i,j = 
0, . . . N — 1. We then determine the direct Lagrangian 
map x(A, t) as the solution of the equation 

d(x,t) = A(x,t) - x . (8) 

First we determine an approximate on-grid solution 
X(A, t) by finding from the inverse map the a point 
nearest to A and its inverse Lagrangian antecedent X. 
We then set x = X + 8x and refine the solution of (jHJ) by 
using a standard Newton method. This requires the cal- 
culation of off-grid values of derivatives, which are again 
obtained from "naive" evaluations of the corresponding 
Fourier series 

P- = ^ifc i d(*.*)e ih -" ) (9) 

OX 3 k 

where d(k, t) are the Fourier coefficients of the displace- 
ment (evaluated in Eulerian coordinates). For each stage 
of the Newton iteration 0(N 4 ) operations are required. 
The number of stages needed to achieve an accuracy e 
consistent with double precision is typically five. If the 
number of output times at which we want to evaluate 
the Lagrangian velocity field is much smaller than the 
resolution N, the displacement-Newton method is much 
faster than particle tracking. 

III. RESULTS 

We have applied the two methods described in the pre- 
vious section to the flow with the initial condition (fTJ). 
The methods give consistent results but the highest res- 
olution (here N — 512) is more easily achieved with the 
displacement-Newton method, which has been used to 
obtain the results reported here. 

The solution of the Euler equation 

du + u-Vu = -Vp, V • u = , (10) 

together with the displacement equation ([7]) was obtained 
by a standard pseudo-spectral method with two-thirds 
dealiasing and a fourth-order Runge-Kutta temporal in- 
tegration. 

Then we applied the displacement-Newton method 
(with five iterations) and e = 10~ 14 . The results were 
checked by computing the Lagrangian vorticity, which 
should be equal to its initial value for 2D Euler flow, and 
was indeed found to be so with an accuracy of 10 -10 . 

In order to locate complex-space singularities for the 
Lagrangian solution, we applied the tracing method Q: 
the Lagrangian solution is represented by its Fourier se- 
ries 

Ma,i)=^^(fc,t)e i,! ' a . (11) 
fc 




FIG. 1: Shell-summed amplitudes of Eulerian and Lagrangian 
velocities at time t = 1.245 in lin-log coordinates. The initial 
velocity is given by ([T}. Inset: time variation (at short times) 
of the width of the analyticity strip in Eulerian coordinates 
(S(t)) and Lagrangian coordinates (<5 L (£)). 

Then the following asymptotic representation is used for 
the shell-summed high wavenumber Fourier amplitudes 

(0) 

I«l(M)I ^C(i)fc a(t) exp[-5 L (i)fc] . (12) 

fc<|fc|<fe+i 

Here <5 L (i) is the width of the Lagrangian analyticity 
strip, that is the distance at time t from the real do- 
main of the nearest (Lagrangian) complex-space singu- 
larity. The same analysis is applied also to the Eulerian 
velocity. 

Figure [1] shows the wavenumber dependence of the 
shell-summed amplitudes for the Eulerian and La- 
grangian velocities at the same time, chosen in such a way 
that there is only a modest range of wavenumbers at high 
k = |fe| at which the rounding errors swamp the (roughly) 
exponential signal. They both exhibit exponential decay 
from which the Eulerian 5(t) and its Lagrangian counter- 
part 5 L (t) are measured. It is seen that the Lagrangian 5 
is significantly smaller than the Eulerian one. Actually, 
S L (t) < 6(t) holds for all times t < 1.245 (the latest time 
analyzed). We also checked that the inequality S h < S 
holds for the flow with the initial condition which has 
less symmetry than (fTJ. 

The exponential decay with the wavenumber k of the 
shell-summed Lagrangian Fourier amplitude is strong ev- 
idence that there are singularities of the Lagrangian ve- 
locity it L (a, t) at a finite distance from the real domain; 
thus it cannot be an entire function. We also obtained nu- 
merical evidence that the Lagrangian map has the same 
locations of complex Lagrangian singularities as the La- 
grangian velocity and that the inverse Lagrangian map 
has the same locations of complex Eulerian singularities 
as the Eulerian velocity. For the very simple Eulerian 
steady flow investigated in Ref. "14,], Lagrangian singu- 
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larities are mapped to Eulerian (complex) infinity. Is 
this also the case for the present flow which has non- 
trivial Eulerian dynamics? Here, the answer appears to 
be "yes". Specifically, let a* be a Lagrangian singular 
location corresponding to time t > 0, say the one clos- 
est to the real domain or one near this position. Does 
<i(a, t) = a x(a, t) tend to infinity as a — > a*? In prin- 
ciple we can find the scaling law of any component of d as 
a — > a*, if we have sufficiently accurate high- resolution 
data for the Fourier transform of d(a) at high wavenum- 
bers. This is explained in Section 4.2 of Ref. [ll|. It 
requires the determination not only of the exponential 
decrement <5 L but of the exponent of the algebraic pref- 
actor in front of the exponential which controls the na- 
ture of the singularity in complex a-space. With a res- 
olution of only 512 2 , such exponents are rather poorly 
determined. It is likely that both components of d(a) 
blow up as s~@ where s is the modulus of a — a+ and the 
exponent f3 is about 3/2 but with and error bar so large 
that a negative value cannot be completely ruled out. 3 
We shall revisit such issues from a theoretical point of 
view in the concluding section. 



IV. CONCLUDING REMARKS 

We have shown that the simple 2D incompressible non- 
steady flow with the initial condition |T]) has complex sin- 
gularities not only in Eulerian but also in Lagrangian co- 
ordinates. The Lagrangian singularities are significantly 
closer to the real domain than the Eulerian ones. A possi- 
ble interpretation of this was given by S. Orszag (private 
communication 2003): in Eulerian coordinates the build 
up of singularities is slowed down by the aforementioned 
phenomenon of depletion, whereas in Lagrangian coordi- 
nates any flow which is non-uniform will keep changing 
non-trivially, even if it is steady in Eulerian coordinates. 
To illustrate this we have shown in Figure the (Eu- 
lerian) Laplacian of the vorticity V 2 a-> in both Eulerian 
and Lagrangian coordinates. The former representation 
displays strongly depleted ribbon-shaped structures, not 
seen in the latter. 

Now we wish to comment on the results concerning the 
analytic structure in Lagrangian coordinates and on pos- 
sible generalizations to other 2D flow with space-periodic 
entire initial data. The most obvious result is that, since 
the vorticity remains unchanged along fluid particle tra- 
jectories in 2D, the Lagrangian vorticity field stays entire 
for all times and thus devoid of any singularities other 
than at complex infinity. The Lagrangian velocity field 



3 We have applied the same method of analysis to the behavior 
of d(x) when approaching a Eulerian singularity at x+. The 
displacement seems again to diverge with and exponent /3 around 
3/2 (implying also the divergence of the Eulerian vorticity) but 
the quality of the scaling is again dubious. 
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FIG. 2: Contours of the Laplacian of the vorticity shown 
for the same flow as in Figure [T] shown at t — 1.245. 
Upper figure: Eulerian coordinates; lower figure: La- 
grangian coordinates. For both figures the contour values 
are 0, ±12.5, ±25, ±50, ±75, ±100, ±125; red: positive, 
green: zero, blue: negative. 



and the Lagrangian map both have complex singulari- 
ties (presumably along one-dimensional complex mani- 
folds) and the numerical evidence is that these are at 
the same locations. Proving this partially can perhaps 
be done by writing the velocity in terms of the vorticity 
using the (periodicity-modified) Biot-Savart integral rep- 
resentation and then making the change of variable from 
Eulerian to Lagrangian coordinates. On the resulting in- 
tegral, using the fact that the initial vorticity is entire, 
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it may be possible to show that if the Lagrangian map 
x(a, t) is analytic for some (complex) a, the same holds 
for the Lagrangian velocity. 

One of the most striking results reported in Section HTTl 
but one for which the evidence is a bit shaky, is that La- 
grangian singularities at time t > correspond to fluid 
particles which at time t escape to infinity. Here are some 
observations which could be useful in proving this. The 
idea is to show that there is a contradiction if at time 
t > a Lagrangian singularity a* at a finite location 
is mapped to a point x± which is not at infinity. In- 
deed, if x± is at finite distance, from the fact that the 
Jacobian of the Lagrangian map is one, it follows that 
x+ must be a singularity of the inverse Lagrangian map 
x i— > a. The Eulerian vorticity can be obtained by com- 
posing the inverse Lagrangian map and the initial (en- 
tire) vorticity. Composing a function singular at x+ with 
one which is entire does not necessarily yield a singular 
function. Perhaps with some extra work it can be proved 
that the Eulerian vorticity is indeed singular at x+. We 
already pointed out in the Introduction that for 2D space- 
periodic initially entire flow there is numerical evidence 
that the vorticity is infinite at (complex) Eulerian singu- 
larities. If this can also be proved, it then follows that 
a* is at infinity and thus we have a contradiction. 

The global picture emerging from all this is (tenta- 
tively) the following: for entire periodic initial data in 2D, 
the solutions of the incompressible Euler equation have 
complex Eulerian singularities corresponding to fluid par- 



ticles initially at infinity and Lagrangian singularities cor- 
responding to fluid particles currently at infinity. In both 
coordinates singularities correspond to some particle es- 
caping to infinity; this mechanism for incompressible flu- 
ids is very different from the one operating for the one 
or multi-dimensional compressible Burgers equation for 
which singularities are mostly associated to the vanish- 
ing of the Jacobian of the the Lagrangian map (see, e.g., 
Ref.[3). 

We cannot at present rule out that the same scenario 
holds in three dimensions but it may not be consistent 
with real blow up. Of course, there are major differences 
in 3D; for example, vorticity is not conserved. However, 
the Lagrangian numerical techniques presented in this 
paper are easily extended to the three-dimensional case. 
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